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Abstract — A  novel  multistatic  array  topology  and  image  recon¬ 
struction  algorithm  for  fast  3D  near  field  microwave  imaging  are 
presented.  Together,  the  techniques  allow  for  hardware  efficient 
realization  of  an  electrically  large  aperture  and  video-rate  image 
reconstruction.  The  array  topology  samples  the  scene  on  a  regular 
grid  of  phase  centers,  using  a  tiling  of  Boundary  Arrays  (BAs). 
Following  a  simple  correction,  the  sampled  data  can  then  be 
processed  with  the  well-known  and  highly  efficient  monostatic 
FFT  imaging  algorithm.  In  this  work,  the  approach  is  described 
and  validated  experimentally  with  the  formation  of  high  quality 
microwave  images.  It  is  further  shown  that  the  scheme  is  more 
than  two  orders  of  magnitude  more  computationally  efficient 
than  the  prevailing  backprojection  method,  and  that  a  cluster  of 
four  COTS  GPUs  can  render  a  3D  image  of  a  human-sized  scene 
in  0.048sec  —  0.101  sec. 

Index  Terms — Microwave  imaging,  multistatic  radar,  Fast 
Fourier  Transform  (FFT). 
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Fig.  1.  1.25m  example  of  multistatic  array  topology  for  fast  near-field 

imaging. 


I.  Introduction 

Near-field  microwave  imaging  is  a  non-ionizing  and  cost  ef¬ 
fective  sensing  modality  for  variety  of  applications,  including 
Non-Destructive  Evaluation  (NDE)  [1],  medical  diagnostics 
[2],  and  detection  of  concealed  weapons  [3]-[5].  In  the  latter 
application,  a  need  exists  for  a  practical  system  that  can  image 
subjects  in  high  foot  traffic  environments,  such  as  mass  transit 
systems,  stadiums,  and  large  public  events.  In  order  to  handle 
a  potentially  constant  stream  of  subjects  in  motion,  the  system 
would  need  to  measure  a  subject  very  quickly  (on  the  order  of 
10’s  of  milliseconds),  and  reconstruct  microwave  images  of  a 
subject  at  video  rate. 

These  system  requirements  present  two  major  challenges. 
The  first  challenge  is  a  cost  effective  realization  of  an  elec¬ 
trically  large  antenna  array.  The  fast  acquisition  requirement 
implies  that  the  system  must  be  fully  electronic  (e.g.,  it 
samples  the  scene  with  no  moving  sensors).  Furthermore,  for 
proper  near-field  illumination  of  a  human  subject,  the  array 
must  be  roughly  the  size  of  the  subject  (l-2m).  For  a  high 
resolution  system  operating  in  the  10’s  of  GHz,  this  implies 
an  aperture  size  in  excess  of  100 A.  A  well-known  paradigm 
that  mitigates  this  challenge  is  multistatic  sampling  [4],  [6], 
[7].  Such  array  topologies  use  transmitters  and  receivers  that 
are  not  co-located.  This  allows  an  array  with  NT  transmit 
elements  and  Nr  receive  elements  to  form  Nr  Nr  spatially 
diverse  samples.  This  is  in  contrast  to  a  monostatic  sampling 
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scheme  (wherein  transmitters  and  receivers  are  co-located) 
which  require  Nr  Nr  transmit-receive  elements  to  achieve  the 
same  sampling. 

The  second  challenge  that  the  required  system  presents  is 
video  rate  image  formation.  The  well-known  backprojection 
algorithm  [4]  can  be  used  with  any  multistatic  configuration, 
but  its  computational  demands  are  extreme.  Fast  Fourier 
Transform  (FFT)  imaging  has  long  been  used  to  efficiently 
construct  images  sampled  with  monostatic  sampling  schemes 
[3],  however,  this  formulation  cannot  be  used  directly  with 
multistatic  sampled  data.  Recently,  a  modified  FFT  imaging 
formulation  for  multistatic  arrays  was  presented  [7].  While 
this  formulation  represents  a  tremendous  improvement  over 
backprojection,  it  still  requires  use  of  multiple  sub-domains 
and  interpolations  to  image  a  human  subject,  adding  substan¬ 
tial  overhead  to  the  FFT  imaging  formulation. 

This  work  presents  a  novel  multistatic  array  sampling 
scheme  that  allows  use  of  the  FFT-based  imaging  algorithm 
after  the  data  undergoes  a  simple  correction.  As  will  be 
discussed,  the  scheme  allows  for  high  quality  image  recon¬ 
struction  of  a  scene  roughly  the  size  of  a  human  subject  at 
video  rate  on  COTS  hardware.  Fig.  1  depicts  the  multistatic 
array  topology.  As  seen,  the  topology  is  a  tiled  arrangement 
of  Boundary  Arrays  (BAs).  The  BA  is  a  well-known  sparse 
array  layout  comprised  of  two  linear  transmit  arrays,  and  two 
linear  receive  arrays  [6] .  A  slightly  different  tiled  arrangement 
of  BAs  was  used  in  [4],  however,  a  key  difference  in  this  work 
is  the  sampling  scheme.  As  will  be  discussed,  the  presented 
design  samples  the  scene  on  a  uniform  phase  center  grid,  with 
no  redundancy.  This  is  critical,  as  it  enables  the  use  of  the 
presented  fast  imaging  technique. 

The  next  section  reviews  the  monostatic  FFT  imaging  for- 
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mulation  upon  which  this  work  is  based.  Section  III  discusses 
the  novel  multistatic  sampling  scheme  and  imaging  algorithm, 
while  experimental  validation  of  the  approach  is  presented 
in  Section  IV.  Processing  load  of  the  imaging  algorithm  is 
discussed  in  Section  V. 

II.  Monostatic  FFT  Imaging 

While  the  presented  imaging  scheme  employs  multistatic 
arrays,  its  explanation  requires  discussion  of  monostatic  FFT 
imaging.  We  consider  the  case  where  a  scene  is  sampled  with 
the  simple  quasi-monostatic  imaging  setup  depicted  in  Fig.  2. 
As  seen,  the  setup  employs  a  raster  scanner  with  two  closely 
spaced  antennas:  one  for  transmission,  one  for  reception.  The 
antennas  are  mated  to  a  Vector  Network  Analyzer  (VNA), 
so  that  the  frequency  response  of  the  scene  is  sampled  on 
a  regularly  spaced  two-dimensional  grid.  Following  collection 
of  all  measurements,  the  image  can  be  reconstructed  with  one 
of  several  techniques. 

Backprojection  [4]  is  the  the  simplest  technique,  where  the 
image  I(r%)  at  voxel  position  is  computed  as: 

NX  Ny 

I{fS)  =EE  s(nx,nv)e^~^M' 

TLX  —  1  Tly  =  \ 

ejk\^  R(nx,ny)-^v\ * 


Fig.  2.  Raster  scanning  setup  for  quasi-mono  static  imaging  and  3D  images 
formed  with  setup  at  18-26.5  GHz. 


It  is  noted  that  s(nx,ny )  is  the  matrix  of  measured  reflections, 
k  is  the  wavenumber,  "TV  is  transmitter  position,  and  r 
is  receiver  position.  As  the  data  is  taken  on  a  2D  spatial 
grid,  nx  and  ny  are  column  and  row  indices  of  the  grid.  If 
reflections  are  captured  at  multiple  frequencies,  then  I(r%) 
is  summed  coherently  across  frequency.  Backprojection  is 
extremely  computationally  expensive,  and  becomes  intractable 
when  the  measurement  domain  is  electrically  large.  FFT 
imaging  [3]  is  a  much  more  efficient  alternative.  Using  this 
approach,  a  2D  image  can  be  formed  with  the  following 
expression: 

I(x,y)  =  IFFT2D[FFT2D[s(nx,ny)}e-jV^^lz °] 

(2) 

where  2:0  is  the  spacing  between  the  measurement  plane  and 
the  image  plane.  Equation  (2)  can,  of  course,  be  evaluated 
at  multiple  frequencies,  and  coherently  summed  within  an 
image  plane.  Spatial  coordinates  (V,  y)  lie  on  the  image  plane, 
and  correspond  to  the  position  of  a  sampled  phase  center 
with  indices  (nx,ny).  For  3D  imaging,  multiple  image  planes 
can  be  evaluated  to  construct  a  volumetric  domain  in  slices. 
Alternatively,  a  formulation  similar  to  (2)  using  3D  FFTs  can 
be  used  [3],  but  this  requires  interpolations  which  may  be 
computationally  expensive. 

Fig.  2  provides  examples  of  3D  images  (formed  in  slices) 
constructed  with  (2).  The  images  were  formed  using  1.25m  x 
1.25m  scan  space,  sampled  on  a  A/^/2-spaced  grid  at  18-26.5 
GHz.  Both  scenes  contain  a  metalized  mannequin  wearing  a 
heavy  winter  coat,  and  a  backpack  concealing  a  can  and  small 
box.  The  images,  which  clearly  show  the  concealed  items,  are 
color-coded  according  to  depth  of  strongest  return,  as  in  [4]. 


The  preceding  example  used  data  collected  with  a  me¬ 
chanically  scanning  setup,  though  an  analogous  electroni¬ 
cally  switching  system  could  be  conceived  for  capture  of 
a  moving  subject.  While  this  would  provide  high  quality 
images  using  an  efficient  reconstruction  technique,  it  would 
require  many  transmit/receive  elements  (roughly  48000  for  this 
1.25m  example),  and  is  likely  not  an  option  for  cost  effective 
array  implementation.  In  the  next  section,  a  hardware  efficient 
multistatic  scheme  is  examined. 

III.  Proposed  Multistatic  Imaging  Scheme 

Multistatic  array  topologies  can  be  used  to  drastically 
reduce  the  number  of  antenna  elements  required  in  a  large 
imaging  aperture,  mitigating  hardware  costs.  Numerous  mul¬ 
tistatic  aperture  sampling  schemes  exist  [4]— [7] .  However, 
with  multistatic  sampling  schemes,  efficient  reconstruction 
remains  a  challenge.  While  the  very  computationally  expensive 
backprojection  method  (1)  can  be  used  for  any  multistatic 
scheme,  the  more  efficient  FFT  method  denoted  by  (2)  cannot 
be  used  directly.  The  multistatic  sampling  scheme  described 
in  this  section  allows  the  FFT-based  method  to  be  used,  after  a 
very  simple  data  correction  is  applied.  As  forthcoming  results 
will  demonstrate,  the  technique  allows  for  video  rate  image 
formation  on  COTS  computing  hardware. 

Fig.  1  depicts  the  array  sampling  scheme.  The  elementary 
unit  of  the  design  is  the  Boundary  Array  (BA)  [6],  a  sparse 
array  topology  first  used  in  ultrasonic  sensing.  This  design 
employs  four  linear  arrays  that  cover  its  perimeter.  In  the 
depicted  example,  the  arrays  on  the  BAs  left  and  right  side 
are  transmitters,  while  the  arrays  on  its  top  and  bottom  are 
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receivers.  The  interior  of  the  BA  is  empty.  Assuming  that 
the  BA  is  sized  such  that  its  side  length  is  short  compared 
to  imaging  range,  each  transmit-receive  element  pair  forms 
an  effective  phase  center,  which  lies  at  the  midpoint  of  the 
two  elements.  If  every  transmit-receive  pair  is  sampled  (e.g., 
via  rapidly  switching  between  the  elements),  a  grid  of  phase 
centers  is  formed.  With  exception  of  the  missing  phase  centers 
in  the  grid’s  middle  row  and  column,  the  grid  is  a  uniformly 
spaced  block.  These  missing  phase  centers  have  a  negligible 
effect  on  image  quality.  In  the  depicted  example,  the  elements 
are  spaced  by  10mm,  or  0.88A  at  26.5  GHz.  The  resultant 
phase  center  grid,  then,  is  spaced  by  0.44 A. 

The  full  array  topology  is  a  tiled  arrangement  of  BAs, 
as  depicted  in  Fig.  1.  Transmit-receive  sampling  occurs  only 
within  a  given  BA  unit  cell:  in  other  words,  a  linear  transmit 
array  never  communicates  with  a  non-adjacent  receive  array. 
The  result  of  this  technique  is  a  scene  that  is  sampled  with  a 
regular  grid  of  phase  centers,  that  covers  the  1.25m  aperture.  It 
is  noted  that  each  phase  center  is  sampled  only  once  under  this 
configuration  (i.e.,  the  topology  samples  with  no  redundancy). 

While  the  sampled  data  set  lies  on  a  regular  phase  center 
grid,  it  still  cannot  be  processed  with  (2)  until  it  undergoes 
a  multistatic-to-monostatic  correction.  We  consider  the  mea¬ 
surements  as  indexed  in  terms  of  the  sampled  2D  grid  of  phase 
centers,  e.g.,  sfax,ny).  Reflections  for  the  missing  phase 
centers  in  the  middle  row  and  column  of  each  BA  are  set  to 
zero.  We  define  a  reference  point  in  the  center  of  the  imaging 
domain,  whose  position  is  given  by  T^re/,  as  depicted  in  Fig. 
3.  The  corrected  data  set  is  given  by: 


sfalx :  Hy)  —  ^fa'x  5 


Rofaxi  Tly) 
Ru  fax  •>  T^y) 


(3) 


Rufai ,  n)  is  the  calculated  reflection  set  that  the  multistatic 
design  would  receive  if  it  probed  a  scene  that  contained  only 
a  point  scatterer  placed  at  faref- 


Rufa'X  77,  )  =  g  jk\^T(nx,ny)  ~^ref  |  g  jk\~^ref  V  R{p>x  iny)\ 

(4) 

R0(m,n )  is  the  calculated  reflection  set  that  would  be 
received  by  a  monostatic  aperture  imaging  the  point  scatterer 
scene,  where  sampling  takes  place  on  the  multistatic  aperture’s 
phase  center  grid  (fac,  see  Fig.  3): 

R0fal,n )  ss=  e-j‘2kfac(nx,ny)-l*’ref\ '  (5) 

The  corrected  data  is  then  zero  padded  to  size  2N  for  FFT 
processing.  In  the  example  depicted  in  Fig.  1,  the  aperture 
supports  a  grid  of  249  x  249  phase  centers.  While  the  data 
could  be  padded  to  256  x  256,  improved  image  quality  was 
obtained  with  a  512  x  512  padding,  and  hence  is  used  in 
subsequent  results.  The  zero  padded  data  is  then  processed 
at  each  frequency  using  (2).  The  process  is  repeated  for  each 
measured  frequency,  and  images  are  summed  over  frequency  at 
each  depth  slice  in  a  3D  imaging  domain.  The  imaging  scheme 
is  summarized  with  the  block  diagram  depicted  in  Fig.  4. 

It  is  noted  that  the  multistatic  to  monostatic  correction  is 
valid  over  a  finite  imaging  domain.  However,  as  illustrated  by 


Fig.  3.  Notional  depiction  of  5x5  tiled  aperture  imaging  a  scene. 


Measured  Fields 


Image 


Monostatic  FFT 
Imaging 


Fig.  4.  Block  diagram  of  fast  multistatic  imaging  algorithm. 


experimental  results  in  the  next  section,  a  single  correction  is 
sufficient  for  a  human-sized  domain.  It  is  also  noted  that  the 
correction  becomes  less  accurate  as  the  size  of  the  BA  tiles  is 
increased. 


IV.  Experimental  Validation 

To  validate  the  sampling  scheme  and  fast  imaging  algo¬ 
rithm,  the  mechanically  scanning  setup  depicted  in  Fig.  5 
was  constructed.  The  system  employs  a  Boundary  Array  (BA) 
emulator:  a  scanner  which  moves  a  pair  of  receiving  antennas 
laterally  and  a  pair  of  transmit  antennas  vertically.  The  four 
antennas  are  mated  to  a  four  port  VNA.  This  allows  for 
sampling  of  the  scene  in  the  same  manner  as  the  BA.  The  BA 
emulator  is  placed  on  an  x-y  scanner,  allowing  for  acquisition 
of  the  tiled  boundary  array  setup  depicted  in  Fig.  1.  It  is  noted 
that,  due  to  the  size  of  the  waveguide  antennas  used  in  the 
setup,  the  closest  transmit-receive  pairs  could  not  be  measured. 
Nevertheless,  high  quality  images  were  still  formed  with  the 
setup. 

Reconstructed  images  acquired  with  the  setup  are  depicted 
in  Fig.  6.  In  one  of  the  two  images,  the  imaged  scene  contains 
a  metalized  mannequin  concealing  a  small  metal  tube  under  a 
heavy  coat,  while  in  the  other,  the  dummy  conceals  a  can  and 
small  box  in  a  backpack.  In  both  cases  the  concealed  objects 
are  clearly  visible.  Both  images  were  formed  with  18-26.5 
GHz  stimulus,  using  160  frequency  points.  The  3D  images 
were  realized  with  21  depth  slices,  spaced  by  0.015m.  Fig.  7 
depicts  the  same  two  scenes,  formed  at  23.6-26.5  GHz,  using 
56  frequency  steps. 
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Fig.  5.  Mechanically  scanning  multistatic  array  emulator  for  validation  of 
sampling  scheme. 


V.  Computational  Load 
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Fig.  6.  Reconstructed  images  at  18-26.5  GHz  (160  frequency  steps)  acquired 
with  setup  depicted  in  Fig.  5.  Formation  time  on  cluster  of  COTS  GPUs  was 
O.lOlsec. 


For  real  time  imaging  capability,  image  reconstruction  must 
occur  at  video  rate.  The  technique  presented  in  the  previous 
section  provides  a  reconstruction  solution  that  is  orders  of 
magnitude  faster  than  backprojection,  and  can  operate  at  video 
rate  using  COTS  computing  hardware.  The  complexity  of  the 
algorithm  was  evaluated  using  computational  workload  values 
provided  in  [8].  The  resulting  load  (in  GFLOPS,  or  billion 
Floating  Point  Operations  per  Second)  is  plotted  in  Fig.  8. 
For  comparison,  the  load  of  the  backprojection  technique  is 
also  plotted.  In  that  analysis,  it  is  assumed  that  only  the 
domain  roughly  occupied  by  the  dummy  (1.9m  x  0.6m  x 
0.3m,  discretized  by  1cm)  is  of  interest.  As  seen,  the  FFT 
technique  is  roughly  two  orders  of  magnitude  less  computa¬ 
tionally  intensive  than  the  backprojection  technique.  For  10  Hz 
operation  with  160  frequency  points,  the  FFT  method’s  load 
is  about  1409  GFLOPS,  while  backprojection  requires  174532 
GFLOPS.  If  56  frequency  points  are  used  at  10  Hz,  the  loads 
of  the  FFT  and  backprojection  methods  are  493  GFLOPS  and 
61148  GFLOPS,  respectively. 

A  COTS  computer  with  four  Nvidia  GeForce  Titan  X 
Graphics  Processing  Units  (GPUs)  was  used  to  reconstruct 
images  using  the  FFT  technique.  For  the  imaging  scenario 
presented  in  Fig.  6  (512  x  512  FFT  size,  160  frequency 
points,  21  depth  slices),  computation  time  was  observed  as 
O.lOlsec.  For  the  56  frequency  point  case  depicted  in  Fig.  7, 
the  computation  time  was  observed  as  0.048 sec.  These  results 
illustrate  that  real  time  image  reconstruction  of  a  human- sized 
domain  is  possible  using  the  presented  technique,  with  COTS 
computing  hardware. 

VI.  Conclusions 
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Fig.  7.  Reconstructed  images  at  23.6-26.5  GHz  (56  frequency  steps)  acquired 
with  setup  depicted  in  Fig.  5.  Formation  time  on  cluster  of  COTS  GPUs  was 
0.048sec. 


highly  efficient  FFT  imaging  technique.  The  method  was 
validated  experimentally,  and  was  shown  to  form  high  quality 
microwave  images.  Computational  complexity  analysis  shows 
that  the  technique  is  roughly  two  orders  of  magnitude  less 
computationally  demanding  than  backprojection.  Further,  the 
scheme  was  shown  to  render  a  3D  image  in  0.048sec-0.101sec 
(depending  on  number  for  frequency  points  used)  on  a  cluster 
of  four  COTS  GPUs. 


A  novel  multistatic  array  topology  and  efficient  FFT-based 
imaging  algorithm  have  been  presented.  The  array  employs 
a  tiling  of  Boundary  Arrays  (BAs)  to  sample  the  scene  on 
a  uniform  grid  of  phase  centers.  After  a  simple  correction 
is  applied,  image  reconstruction  can  be  performed  using  the 
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